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Abstract 

We study the asymptotic entanglement of three identical qubits under the 
action of a Markovian open system dynamics that does not distinguish them. 
We show that by adding a completely depolarized qubit to a special class of two 
qubit states, by letting them reach the asymptotic state and by finally eliminating 
the added qubit, can provide more entanglement than by direct immersion of the 
two qubits within the same environment. 

1 Introduction 

A quantum system S interacting weakly with its environment E can be treated as 
an open quantum system: standard techniques can be apphed to obtain a master 
equation involving the degrees of freedom of S, only. On a long time-scale determined 
by the weakness of the coupling to the environment, the reduced dynamics consists of a 
semigroup of trace-preserving completely positive maps on the states of S. These maps 
are fully consistent transformations of the system S states incorporating the dissipative 
and noisy effects due to the environment -E [H [2] • 

An irreversible reduced dynamics typically transforms pure states into mixed states, 
thus spoiling fundamental quantum resources like entanglement p]. However, the 
Markovian reduced dynamics resulting from suitably engineered environments may 
entangle an initial separable state of a bipartite quantum system immersed within it; 
such an entanglement can even persist asymptotically IH El El 13 [8l [9l [101 [HI 112] ■ 

The properties of the asymptotic states of quantum dynamical semigroups can be 
studied by looking at the structure of the generator of the reduced dynamics [IH [HI 
[151 [m [m [m [19]: some concrete applications can be found in [9l [121 [201 [211 [22] . 



Most of the results present in the hterature concern pairs of open qubits; in the 
following, we will instead consider three qubits weakly interacting with a common en- 
vironment such that the resulting master equation affects all possible pairs of qubits in 
the same way. In this case, the manifold of stationary states can partially be analyti- 
cally characterized as well as the asymptotics of a particular class of initial states. We 
shall use the analytically obtained results to construct the following protocol: given 
two qubits immersed in the just depicted environment, 1) we add a third completely 
depolarized qubit, 2) leave the three of them irreversibly evolve until they reach the 
stationary state, 3) trace away the ancilla. As an indication of the surprises set in store 
by higher dimensional entanglement with respect to the two qubit case, we show that, 
for certain initial two-qubit states, the entanglement obtainable by such a procedure 
is larger than that achievable by simply letting them reach the stationary state within 
the bath. 



2 Three open qubit model 

In the following we shall study an open quantum system of three qubits in weak inter- 
action with their environment E. We shall denote by M the algebra of 8 x 8 matrices 
X G Ms{C) and by p G M the positive matrices of trace 1 that describe the states of 
5* and by S{M) the convex set of all states. For sake of simplicity, we shall sometimes 
write p{x) for the expectation values Tr(px); further, by x^^\ y^'^\ z^^\ respectively 
Pi^^,p^\pf\ we will denote the local one-qubit observables x®!®!, l^ycg)!, 
1 ® 1 ® 2;, respectively the one-qubit states pi®!®!, 1®P2®1 and 1 (g) 1 (g) P3. 

In the weak-coupling limit, the states of an open quantum system evolve in time 
according to a master equation of the form dtPt = L[pt], where the generator L incor- 
porates the dissipative and noisy effects due to the environment; the solutions form a 
semigroup of completely positive maps 7* = exp(tL), t > [H [21 [9] . In the present 
case of three qubits, we shall concretely study the following time-evolution equation, 

K"" = [ E Y -3°' . + E (-."" - 5 {-f'-.'°' . ft}) = . (1) 

a=l a, 6=1 

where a^""^ is the i-th Pauli matrix of the a-th qubit, and the coefficients C^^j''^ form a 
a 9 X 9 positive matrix, the so-called Kossakowski matrix: 

/ (7(11) (7(12) (7(i3)\ 

To the semigroup of completely positive, trace-preserving maps = exp(tL) on the 
states of S, there corresponds the semigroup of completely positive, identity-preserving 
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maps 7j : M I— )■ M, on the matrices x G M, generated by: 



dtXt 



6 
a=l 



'a (a) 



a, 6=1 



(ab) 



(a) 



(fe) (a) 



L[xt] . (3) 



The analytic solution of ([T]) can in general be addressed only numerically; we shall 
instead consider a simpler class of master equations amenable to a partially analytic 
study. Notice that, by taking the trace of the generator L in ([1]) with respect to any 
single qubit, one gets the generator of a master equation relative to the other two 
qubits; thus, by choosing ui = ijJ2 = 0U3 = ou and 

^ (j{22) ^ (j{33) ^ ^ . (j{12) ^ (j{2l) ^ (j{13) ^ (j{31) ^ (j{23) ^ (j{32) ^ ^ ^ 

one obtains a highly symmetric generator with Kossakowski matrix 



K 




(4) 



such that any pair of qubits is affected in the same way by the presence of the environ- 
ment. As apparent from (jll), the matrix A governs the dissipative reduced dynamics 
of each one of the qubits, while B rules the dissipative statistical coupling of pairs of 
different qubits. However, in the following, we shall further restrict the case to master 
equations where one and two-qubit terms are the same and choose 



A = B 




c > , < |6| < a . 



(5) 



Then, the Kossakowski matrix ([2]) reads K = 3 where P projects onto the vector 

(1,1,1)/a/3 and A > 0. The resulting master equation for the states (Schrodinger 
time-evolution) and its dual for the system operators (Heisenberg time-evolution) can 
conveniently be recast as 



dipt 
dtXt 



S3, pt 

S3, Xt 



3 ^ 

^ij (^Sj Xt Si — — ^^Sj Si , Xt^^ = L[xt] 



in terms of global spin operators 

3 



Si 



2,3 



a=l 



(a) 
2,3 



(1) 

1,2,3 



+ a 



(2) 
1,2,3 



+ a 



(3) 
1,2,3 



(6) 
(7) 



(8) 
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Remark 1. For three qubits (n = 3), a master equation of the form 

3 



dipt 



u 
-I — 
2 



(a) (a) 
ij \ Pt CTj 



a=l i,j=l 



(a) (a) 



J * ' 



, (b) (a) 

- i a) 'a\ , Pt 



}) 



aj^b=l i,j=l 

may have direct experimental imphcation in certain reahzations of the driven cavity 
array proposed in [23]. For n = 2, the above equation have been derived in a physical 
scenario where the qubits are at a distance from each and immersed in a scalar Bose 
field in thermal equilibrium [20], while a master equation of the form 

3 3 



dtpt 



''2 



S3 , Pt 



a, 6=1 i,j=l 



(b) (a) 



corresponds to two qubit immersed in an environment described by a thermal, scalar 
Bose field when the spatial distance among the qubits is negligible [9]. 



3 Asymptotic States 

We shall start by briefly reviewing some available results about the stationary states 
of quantum dynamical semigroups ||T3| [HI [15] and about the tendency to equilibrium 
of open quantum systems. 

Let iS^ = {p G iS(M) : ■jt[p] = p Vt > 0} denote the set of stationary states of a 
semigroup of trace-preserving, completely positive maps ■jt = exp(tL), generated by 
the master equation dtPt = L[pt], and by = {x E M : 7t[x] = x Vt > 0} the set of 
operators {n x n matrices) invariant under the identity-preserving maps ^t = exp(tL) 
generated by the dual time-evolution equation dtXt = L[xt]. It is convenient to cast 
the latter equation in diagonal form: 

3 , 

+ 5^(\//a;tK. . (9) 

i=l 

From [15] one knows that, if a full-rank stationary state poo exists, then 

1. the subset of constant matrices = {x E M : 7t[x] = x} is a *-subalgebra of 
M, that is ^tix"^] = x"^ and also 7j[x'''x] = x^x for all t > 0. 

2. the time-average 

1 

E[x]= hm - / dt%[x] (10) 

defines a conditional expectation from M onto M^, that is a completely positive 
unital map such that 

E[l] = l, t\yixy2]=yit[x]y2 Vt/i,2 G , Vx G M . (11) 



dtXt 



if, Xt 
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The conditional expectation E : M k-)- has a dual map defined by 

Tt(^pE[x]^ =Tt(e[p]x^ , \f p e S{M) , X e M . (12) 

This is a completely positive, trace-preserving linear map on the state-space S{M) 
such that E[p] is a stationary state and E[p] = p if p is a stationary state. 

We are interested in establishing whether, given any initial state p, it goes into an 
asymptotic state poo according to 

Poo= lim 74p]=E[p] . (13) 

i— >-+oo 

A sufficient condition can be obtained as follows: consider the subset C M of 
X G M such that 

I/[a;^a;] — x — a;''"I/[x] =0 . 
From ([9]) one immediately derives that C and also that x G if and only if 

3 

5^([x, y,]) + [x, V^ = 0^[x, Vi\ = VF, . 

i=l 

Thus, the subset consists of x G M which commute with all operators Vi, namely 
= {Vi}' where {V^}' denotes the so-called commutant of the set {Vj}. The com- 
mutant is a subalgebra which need not coincide with the time-invariant *-subalgebra 
My. This is however the case if the operators [15] commuting with all Vi also commute 
with their adjoints and with the Hamiltonian H. Indeed, if {Vi}' = {Vi,V^^ , H}' , then 
M^CD^C as ([9]) yields {V, 1//, H}' C M^. 

Moreover, the equality {Vi}' = {Vi,Vi , H}' is also sufficient [15] to guarantee that 
lim 7t[x] = E[x], for all x G M, whence f|T3l) follows by duality. 

In general, that is for any number of qubits, the time-evolution equation ([7]) can be 
written as in by diagonalizing the 2x2 matrix in the upper left corner of A in ([S]); 
concretely, in terms of the spin operators in ([H]), 

Vl, = V2(^ ^^^^ , V = V~cSs. (14) 

If \b\ < a and c > 0, {Vi}' = {V, V^ , H}' = {Si}' so that, according to the above dis- 
cussion, it follows that = {Si}'. Therefore, in the concrete cases we are considering, 
the time-invariant operators coincide with those commuting with all global spin oper- 
ators in ([8]). In order to establish the asymptotic convergence to stationary states 
as in f[T51) . we need seek full rank stationary states: we shall do this in the following 
for 1, 2 and 3 qubits. 



5 



3.1 One qubit 

For the case of one qubit, a full rank stationary state of the master equation 

3 ^ 
+ ~ 2 {^^ ' ^*}) ' ^'^^^ 

with A = [Aij] as in (E]), can be found by considering the corresponding time-evolution 
equation of the Bloch vector ft in Pt = -(1 + rt ■ (?); that is, f = —2{Cft — z), where 

/a + c -uj/2 0\ /0\ 
C = I uj/2 a + c , i'= . (16) 
\ 2a/ \2bJ 

Setting ft = 0, one finds foo = Cr^z = {0,0, b/a) and a unique full-rank (|6| < a) 
stationary state 

P*oo = 2(^ + ^00(^3) , ^°°=a- ^^^^ 

Remark 2. Consider two qubits (a = 1,2 in ([8])); one explicitly verifies (see also [9j) 
that p®^ = ® p;!;^ is a full-rank stationary state for (E]): L[p'^] = 0. 

The generator of ([6]) can be extended to the case of n qubits by extending to n the 
summation index of single qubit Pauli matrices in (j8]); further, it can conveniently be 
recast as L = ^"^,=1 Lab where the sum is over generators ([6]) involving only the ath 
and 6th qubit. Let p®" = p*^ p*^- ■ ■ p*^ ; then, 

n times 

L^''^[pT] ■■= {Lu + L,, + + L21) [pT] = l^'IpT] ® p^(""^) = , 

where L^^ is the generator in ([6]) for two qubits and p^^ is a two qubit stationary state. 
This result clearly holds for all pairs (ab), that is L("''^[p®"] = 0, whence L[p^] = 
and p^" is an n-qubit full-rank stationary state. 

3.2 Two qubits 

As previously observed, {Vi}' = {Vi, , H}' = {Si}' independently of the number of 
qubits. In order to find the commutant {Si}' for the case of two qubits, we use the 
Pauli matrices and write 

3 2 3 
i=l a=l i,j=^ 

Then, by imposing that [x , Sp] = for all p = 1,2,3, {Vi}' amounts to being the 
linear span of the identity matrix 1 and of the symmetric sum T = Yl^=i'^i ® 



dtpt = -i 



0"3, pt 
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It follows that = {Si}' is a commutative algebra; it coincides with its center, 
= Z = {Si}' n{Si}" = MjHM!^ and is generated by the two orthogonal projections 

P=i(l-T), Q = i_P=i(3 + T) , (18) 
where the first one is 1-dimensional and projects onto the two-qubit singlet state 

|^) = -^(|0)®|1) - |1)®|0)) , (19) 

with (T3IO) = |0) and a^ll) = -|1). 

From Remark 2., p®^ = (g) pj!j^ is a full-rank stationary state; then, f fT3l) ensures 
that the asymptotic state poo corresponding to an initial p is obtained as IE[p], by means 
of ( fT2l) . In order to construct it, we first construct the conditional expectation E onto 
the sub-algebra of constant matrices; E must be such that E[x] = A(x) P + p(x) Q. 
From the properties (ITTi) of the conditional expectation, 

E[PxP] = X{x) P , E[QxQ] = p(x) Q , 

where, with p(x) := Tr(px), 

_ Tr(Ppf Px) _ Tr(gpf gx) 

Then, from (IT^ one gets 

pf^=^^^', roo = -, (20) 
4 a 

so that, given any initial state p, its asymptotic state poo is given by (compare with [H]) 
Poo = E[p] = 



The entanglement content of any two-qubit state p is quantified by the concurrence 
C{p) [21]: consider the complex conjugate matrix p*, construct p = cr2®cr2 p* cr2^cr2 and 
compute the (positive) eigenvalues Xj of pp. Then, C(p) = max{0, Ai — A2 — A3 — A4}. 
For all asymptotic states poo in fl2T|) . one easily calculates 

C^(Poo) = ^(^^7^ max {0, 2 |4p(P) - (1 - r^) I - 2(1 - p(P))(l - r^)} (22) 

In [9] the entanglement capability of the environment has been studied by comparing 
the concurrence of certain initial states with that of their asymptotes; in the following, 
we shall focus upon the following one-parameter family of initial conditions 

p{a) = al + (l-4a)P, 0<a<l/3. (23) 



4p(P) 
1 — tH-, 
4(1 



4p(g) 

3 + r2 



QpZ'Q 



3 + r2 



Poo ' 



4p{P) - 1 + r 
3 + r2 



2 

00 p 



(21) 
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One easily finds that C{p{a)) = max{0, 1 — 6a}. Furthermore, if 



< « < «(r^) = f '\ , (24) 
6(3 - r^) 



where a^r^o) is an increasing function of Too: 1/6 < a(roo) < 1/3, the corresponding 
asymptotic states obtained, according to f|T3|) . as 



= E[p(a)] = + P (25) 

-J "T ' oo "J "T ' OO 



have concurrence 

C{Poo{a)) = l-3a^^ > 0. (26) 

^ "J 'oo 

Otherwise, namely for a{roa) < «, Poo(«) is separable. One can then conclude: 

1. both p{a) and Poo(tt) are separable if 

I < «(roo) < « < ^ . (27) 

6 6 

2. p{a) is separable and Poo(«) is entangled if 

^ < a < a(roo) ■ (28) 
6 

3. Since a{roo) > 1/6, it follows that, when < a < 1/6, the initial state p(a) is 
entangled as well as poo(a); the entanglement difference 

A(«) := C(poo(«)) - C{p{a)) = 9 a - ]- (29) 

"J + ^oo 



becomes positive (entanglement gain) if 



« > «*(^oo) = J.^'h, ' (30) 
18(1 + r^) 



where a*(roo) is a monotonically decreasing function of Too: 1/6 > a*(roo) > 1/9. 



3.3 Three qubits 

As for one and two qubits, in order to fully characterize the set of asymptotic states, 
one needs the conditional expectation ffTOj) : differently from two qubits, in the case 
of three qubits its complete expression is still escaping us. Indeed, the commutant 
is not commutative and cannot coincide with its center, ^ Z (see Appendix 
A); neither does coincide with the commutant of its center (M^ ^ Z'\ which is 
the other case where one would immediately know how to construct the conditional 
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expectation [131 US] • What can be analytically constructed is at least the action of E 
on certain subsets of initial states. 

In appendix A, it is showed that the commutant set {Si}' = is the linear span 
of the the 3x3 identity matrix and of the following operators 



= J2 -^-f ^ a < 6 = 2, 3 ; S = a^^afa^ . (31) 

i,j,k=l 



1=1 



Further, the center Z surely contains the operator 



T= J2 3'^'' = 5(12) ^ ^(23) ^ ^(13) _ 



a<b=2 



Also, the operators 



piat) ^ ^ ^ ^ M 

4 ^ 



(32) 



(33) 



are projections such that [p("^) , Si] = for all i = 1,2,3. Using ([9]), given any 
p G S{M), the states 

p{ab) p p(ab) 



Jab) 



are such that 



l^piab)^ = pi-b) ^[p] p{ab) ^ ^^[^(a6)| ^ p(afe) p(af,) _ 



(34) 



(35) 



Moreover, as P^"''^ = \'^ab){'^ab\ Ic projects onto the singlet vector state {"^ab) of the 
qubits a and b, then 

lt[Pab]=P^^'^P^^ , (36) 

(c) 

where is a state of the qubit c. 

Proposition 1 The state p['^^ evolves in time according to the master equation ([6]) 
for one qubit and E[p('^^)] = P("''V^\ where p^^ = p*^ in (fT7|). 

(c) 

Proof: The time-evolution of pi is obtained by tracing over the qubits a and b the 
expression (|36|) multiplied by p("^); by using fl35|) one gets: 

5,pl^)=Tr„,(p('^^)L[p(«'')pl^) 
By splitting the generator as L = ^^-^ Lpg, one gets 



Lii 

+ P^^'^L^pP]. 
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The first contribution vanishes for it consists of the generator of the master equation ([6]) 
for two qubits acting on the projection onto the singlet state; from (^^, this state is 
stationary and the statement follows. 

Since P^"''^ e {Si}', the trace over the qubits a and b of the second contribution 
multiplied by P^"'') reads 



This piece vanishes, too; indeed, all Kraus operators contribute with terms of the 
form Tr(p("^^ crl"'') or Tr(P^"'') '^f'^) which are both zero as the partial trace of P^'^'') is 
proportional to the 2x2 identity matrix and the Pauli matrices are traceless. Therefore, 
dtPt = LcclPt^'^] whence the result follows from the fact that Lcc is the generator in ([6]) 
for a single qubit which has as full rank stationary state. 

!_ _ 5" (aft) 

One can now fix the action on the projectors P^"-''^ = ^ and 

P=- J2 ^^"'^ (37) 

a<b=2 

of the dual map E introduced in ( !T2l) which, according to (fT3ll . associates to any initial 
condition the asymptotic states towards which it tends when t — )■ +oo. 



A 3 

4 "^-.,6) (c) 

r oo • 



Corollary 1 E[p("'')] = 2P('^''V^^ and E[P] = - ^ P^"' 

a<b=2 

Proof: Set p = p("^) in then, E[P^''^'^] = lim 7t[p(°''^] = 2 P^"''^ p'^^ . The second 

t^-+oo 

relation follows by the linearity of E. 

p{ab) 

Remark 3. Notice that while P("^) G and thus E[p("^)] = p("^), is not an 



invariant state: E 



p{ab) 



p{ab) 



The last necessary tool for the applications to be discussed in the next section is the 
action of the map E on the projection Q = \ — P ^ M^. 



Proposition 2 E[Q] = -^(f' - E ^^"'^Z'- 



Proof: Since Q G M^, the properties (jlTi) and the algebraic relations (15511 applied to 

E[x] = \{x) 1 + Y.l<b=2 ^-h{x) S^^'^ + p{x) S,xeM, give 

3 

E[QxQ] = QE[x]Q = P{x)Q , f3{x) = X{x) + ^ Xab{x) . 

a,b=l 
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Using the time- invariant state p'^, IE[p^^] = p'^, one obtains 

Tt(QpZ'Qx) = Tr(pf QxQ) =Tr(E[p®3]QxQ) =Tr(pf E[gxQ] 
= Trfpf E[Qxg]) =/3(x)Tr(pf g) . 



Tr(gp^3ga:) 

This gives /3(x) = — ; on the other hand, for all x G M, 

Tr (pgf g) 

Tr(a;E[g]) =Tr(E[gxg]) = ;^r|P^Tr(xgpf g 
Then, the result follows using that (see ( 120|) ) 

a<6=2 a<b=2 

In the next section, we study the following protocol: 

• add a third completely depolarized qubit to a two qubit initial state p(a) as 
in (El; 



• let the resulting three qubit state reach equilibirum under the time-evolution 
governed by ([6]); 

• eliminate from the asymptotic state the added third qubit. 
We show that the resulting two qubit state 

1. can be entangled when the asymptotic state reached by the two qubits evolving 
alone would not, that is when a > a{roo); 

2. can be more entangled than the initial state, when the asymptotic state of the 
two qubits evolving alone would not, namely when a < a*(roo); 

3. the entanglement gain can be larger than A (a) > when a > a*(roo)- 



4 Applications 

We now apply the previous results to the study of the asymptotic entanglement proper- 
ties of a class of three-qubit states obtained by appending to the two-qubit states f l23|) a 
third qubit in the completely depolarized state; we shall thus focus onto initial density 
matrices of the form 

pi23(«) = p(«) ® 1 = ^ 1 + 1^1^ p(i2) ^ < a < 1/3 , (38) 
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where, according to the notation of the previous section, P (8> 1 = P*-^^^ 

The corresponding asymptotic states are given by the map E : S{M) i— )■ S{M) whose 
action is given by Corollary 1 and Proposition 2; indeed, writing 1 = P + Q, 

m = yt^pT + ^^^^ E (39) 



a<6=2 



+ (l-4a)p(i2)pS^ (40) 

According to the last step of the protocol described at the end of the previous section, 
we trace the asymptotic states p^^(a) with respect to the appended qubit: 

19/ N / 19^/ XN 4a 4ar^ + 3(1 - 4a)(l + ) 



1 + 3(l+r2 



OOJ 



2 Of / \ 

+ TTTT— ^ 1 ® Poo + Poo ® 1 , (41) 

3(1 + ^^4) ^ ^ 

where P projects onto the two-qubit singlet state. The concurrence C{p^a?\a)) of this 
two-qubit state can be computed and compared with that of the asymptotic state p°°(a;) 
in (l30l) : though easy to calculate, the expression of the concurrence is not particularly 
inspiring and can be found in Appendix B, equation (1581) . The goal is to see whether 
the addition and final discarding of the added qubit may increase the asymptotic 
entanglement of Poo(a) in fl25l) . 

We start by considering the case of separable two-qubit state p(a) that cannot get 
asymptotically entangled by the action of the master equation ([6]). According to (!27|) . 
this occurs for 1/6 < a(roo) < « < 1/3. 

Consider a third qubit prepared in the totally depolarized state and appended to 
the qubits 1 and 2 prepared in a state p{a) with a in the above range. According to 
Appendix B, by tracing the asymptotic 3-qubit state p^^(a) over the appended qubit, 
the qubits 1 and 2 are entangled, that is their concurrence C{p^'^{a)) > 0, if either 
< a < a;+(roo) where 



, , 3(1 + r 



4{3 + 2ri^) + 2./6{^ 



or a_(roo) < « < - where 
3 



. X 3(1 + r 



4(3 + 2r^)-2v/5(^ 

with (5(roo) = (1 - r^)((3 + 2r^)2 - 9r^) > 0. One checks that a+{r^) < a{r^y, 
therefore, the first condition is incompatible with (12T 
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We shall then set so that a.^r^) < a < - and let 0.980965 = r* < < 1 

as calculated in Appendix B, equation (!64l) . If a_(roo) < a(?"oo), all initial states 
p{a) with a > a(roo) correspond to asymptotic states Poo(tt) which are separable, 
but to reduced asymptotic states pl^{a) in (HTl) that are entangled. The same occurs 
for a_(roo) > Q!('"oo) for initial states p{a) with a > a;_(roo). Therefore, there are 
separable states p{a) which do not get asymptotically entangled by direct immersion 
in the environment described by IQ, but do get entangled if a third depolarized qubit 
is appended to them and then eliminated after reaching stationarity. 

This phenomenon is numerically confirmed in the following figure where the concur- 
rence of p^{a) is plotted for < Too < 1 and 1/5 < a{roo) < a < 1/3. 



0.04, 
0.03- 
0.02- 




Figure 1: C(p^(q;)) when C(poo(a)) = C{p{a)) = 0, r = Too, a = a 



Remark 4. Admittedly, the range of favorable values of the environment dependent 
parameter Too is not so large, as well as the range of separable two-qubit states p{a) 
that can get entangled by means of the protocol and not by direct immersion in the 
environment. However, the fact that such a possibility exists is an indication of what 
might be achievable if one could completely characterize the whole manifold of sta- 
tionary three-qubit states. Also, instead of tracing away the third qubit, one could 
perform a less mixing operation on it in such a way that some more entanglement be 
localized on the remaining two qubits: preliminary results confirm this possibility, but, 
unfortunately, not to a sufficiently significative extent. 

Luckily, concerning the second two points listed at the end of Section 3.3, addition of 
a third completely depolarized qubit and its elimination after reaching the stationary 
regime, allows for a more substantial improvement on the entanglement that can be 
gained asymptotically. Let us consider the difference 

Ai(«):=C(p^(«)) -C{p{a)) (42) 

in the range a < a*{roo)- For these values of a, no entanglement gain can be achieved 
by letting the two open qubits evolve towards their stationary state; that is, 

A{a)=CipM)-Cipia))<0 . 
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However, by adding a completely depolarized qubit and eliminating it after reaching 
the stationary state, one may get Ai{a) > as showed in Figure 2 which exhibits the 
range of parameters r^o (depending on the environment) and a (labeling the initial 
state) for which this is possible. 




Figure 2: Ai{a) vs r — r^, < a — a < Q;*(roo) < 1/6 

Next, consider the difference 

A^ia) := C(p^(a)) - C(p^(a)) (43) 

in the range a > a*{r^) where two qubits present an entanglement gain, A(a) > 0. 
Such an entanglement gain may be increased by adding a third depolarized qubit as 
showed in Figure 3. 




1 0.1 



Figure 3: A2(a) vs r = Too, «*(?'oo) < a = a < 1/3 



5 Conclusions 

We have studied the asymptotic states of a Lindblad master equation describing the 
reduced dynamics of three qubits weakly coupled to an environment that affects in 
the same way any pair of qubits. By applying standard algebraic techniques, we could 
control the asymptotic states of a particular family of initial three qubit states of 
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which one is completely depolarized. We showed that, after eliminating the latter 
from the asymptotic state, the remaining two qubits may show more entanglement 
than the asymptotic two-qubit state achievable by direct immersion within such an 
environment. This phenomenon can be regarded as an asymptotic manifestation of the 
richer structure of irreversible entanglement generation in higher dimensional discrete 



systems that was observed at short times in [10 



Appendix A 

Given the operators Si = Yll=i "^ = 15 2, 3, the commutant set = {Si}' is found 
by expanding a generic x G M by means of tensor products of the Pauli matrices: 



^ - \ 1 ^ \^ X^'^K'-"^ J- Sr \(2) ^W^W I \^ \(3)^(l)^(2) 



(3) 
k ' 



(44) 



a=lA=l 



a<b=2;i,j=l 



i,j,k=l 



and then imposing [x , Si] = for all i = 1,2, 3. By using the Pauli algebraic relations 
one finds the following equalities 

3 

J2^ae^ipi = ^ Vz,p=l,2,3 (45) 

1=1 
3 

E (^Sfe£ + X%Jeep^) =0 V a < 6 = 2, 3 ; z, j, p = 1, 2, 3 (46) 
3 

[^ifi ^tpk + ^tpj + >Hjk ^ipi) = Vz,j,A;,p=l,2,3, (47) 

1=1 

whence A^^'* = for all a,i = 1, 2, 3, A^^^^^ = A^^\^- for all a < 6 = 2, 3 and i,j = 1, 2, 3, 

while Xl^l = Xsijk- It thus follows that the commutant set is {Si}' = {1, S^'^^\ S}, 
a,h = 1, 2, 3, namely the linear span of 1 and 

^(-) = ^ ,(^) ,f ) a<b = 2,3; S = e.,k 4^ of . 



(4J 



i=l 



i,j,k=l 



Unlike for two qubits, the commutant set is not commutative; indeed, with a, b, c 
different indices, 



S =42(^^('"=) , a<b 



(49) 
(50) 



whence 



T = J2 S'''' = 5(^2) + S^^^^ + S^'^^ [T , S^'''^] = [T,S] = , (51) 



a<b=2 
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so that T belongs to the center Z = {Si}' fl {Si}" = 
relations are as follows 

(^H))2 = 3_2^('^^) , a,6 = 1,2,3; 

From the first relations it follows that 

1 - S^"'''^ 



My n M'. Other useful algebraic 



S^ = 2(3 - T) 



p{ab) 
P 



9 ^ 

2 y 

3 ^ 

a<b=2 



eM^ = {Si}' = 



p(ab) ^ _ ^rp^ ^ ^ 



(52) 

(53) 
(54) 



are two-dimensional, respectively four-dimensional projections. In particular, P^"-^^ is 
the tensor product of the projection onto the singlet state of the qubits a and b with 
the identity matrix for the qubit c. Furthermore, the projection Q = 1 — P E fl 
fulfils 

gS^'^^) = Q ya<b; QS = 0. (55) 

Other projections commuting with M^, that is in the commutant are thus all 
sub-projections q < Q for which qQ = Q = Qq, whence 

g^C'^^') = gQ5("^) = gQ = g , qS = qQS = 0. 

However, unless q = Q, these projections q cannot belong also to M^; this is proved by 
writing g = Al -|- ^a<;,=2 -^af' and by imposing the previous conditions. It 
thus follows that, for the three-qubit case discussed in this work, neither Z = as 
this would imply commutative, or Z' = as this would imply C M^. These 
are the two conditions for which a conditional expectation onto could easily be 
explicitly written 



Appendix B 

The explicit form of the states p^(a) in fj4T]) is 



where 



u 



2 

CXD 



m 


is 














o\ 





y 


—u 








—u 


y 





\o 











(56) 



a 



-(l±roo)(3(l± 
3(l + r^) -4a (3 + 2r 



2ri 



3(l+r^)-2a(3 + 5r2 



oo/ 



6 



(57) 
(58) 
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The concurrence of such a state is 

-2aV(l-r^)(9 + 9r^ + 4rl)| . (59) 
More exphcitly, set 6{r^) = (1 - r^) (^(3 + 2rl^f - 9r^j > 0; then, 

= ^ ^^^^ ^-1 (60) 



3" ' 4(3 + 2ry-2y5(^ 

Cip-ia)) = 1 ^-^^^ ^ (61) 

if < a < a+{r^) - ^ ' co) 



4{3 + 2rl) + 2,/S^ 
The lower bound Q;_(roo) is a decreasing function, 

A = «_(1) < a_(roo) < ^ = a_(0) , (62) 

while the upper bound a^{r oo) monotonically increases, 

^ = a+(0) < a+(r,,) < ^ = «+(l) . (63) 

While a_|_(roo) is always in the permitted range < a < -, it turns out that 

(y-{roo) < ^ if 0.980965 = r* < Too < 1 , (64) 
where r* is such that Q;_(r*) = 1/3. 
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